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Abstract 

The electroweak corrections of order ©(ocm'Tij j/m^) and 0{aauim^ fj/m^) to Higgs boson associated 
production with top quark pair are calculated at e'^e~ colliders in the standard model (SM), the two-Higgs- 
doublet model (2HDM) and the minimal supersymmetric standard model (MSSM). These corrections are 
a few percent in general, and in the MSSM they can be over ten percent for favorable parameter values 
allowing by current precise experiments. The total cross sections including the electroweak corrections get 
their maximal near ^ = 700GeV, and can reach 2.8 fb, 2.7 fb and 2.5 fb in the SM, the 2HDM and the 
MSSM, respectively. 
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1 Introduction 



Higgs mechanism [1] is a milestone to the establishment of Standard Model [2] (SM), where the SU{2)l Higgs 
doublet generates electroweak symmetry breaking, gives masses to the weak gauge bosons W^, Z through the 
kinematic terms of Higgs fields, and gives masses to fcrmions through Yukawa coupling of Higgs with fermions. 
In the SM with one SU (2)^ Higgs doublet, there is only one physical scalar Higgs boson left after the electroweak 
symmetry breaking, the CP even Higgs boson /i**™, while in the two-Higgs-doublet model (2HDM) with two 
SU{2)l Higgs doublets, there arc three neutral and two charged Higgs bosons /?,, H, A. and H^, of which h and 
H are CP even and A is CP odd. Due to a fundamental scalar /t"™, there exists "hierarchy problem" in the SM, 
where the Higgs self-energy is quadratic divergent from loop of Dirac fermions, which pulls the Higgs mass to a 
high cut-off scale [3] . Supersymmetry is introduced to solve this problem, through the cancellation of quadratic 
divergence in the loops of Higgs self-energy, in which 2HDM of usual type II have to be introduced [4] , where 
one Higgs doublet couples to weak isospin I = \ fermions, and the other Higgs doublet couples to / = — ^ 
fermions. 

The Search for the Higgs boson is one of major objectives of future high-energy collider, after the discovery 
of top quark [5]. Great effort has been devoted to this direction to study different properties of the Higgs boson. 
The direct search of Higgs boson at LEP H through the process e+e" Zh^™ has set the lower experimental 
limit of 114.1GeV on the Standard Model CP even Higgs boson mass [6]. Indirect constraint on the mass of 
Higgs boson from the W mass measured at Tevatron and LEP II and top quark mass measured at Tevatron 
is below 195GeV at 95% confidence level [7], which signifies a light Higgs boson. The latest lower limit on the 
lightest CP-even Higgs boson mass in MSSM is ruh > Ql.OGeV [8, 9], where production process e+e~ — > Zh is 
suppressed with small sin^(/3 — a). When sin^(/3 — a) is not suppressed, the limit on the lightest CP-even Higgs 
boson mass from above process is about lI4GeV [9] in the MSSM. While in the MSSM, because supersymmetry 
relates Higgs quadratic coupling to the gauge coupling, there exists an upper limit on the lightest CP-even 
Higgs boson mass m/j < ISOGeV [10] when including radiative corrections. 

With a heavy top quark, tth associated production also can be helpful for both the discovery of a Higgs 
boson in the intermediate mass range and the measurement of large top quark Yukawa coupling at the near 
future colliders. There are the possibilities at the Tevatron Run II to search the Higgs boson through qq tih 
production subprocess, followed hy h ^ bb, for < 140GeV [11], though the statistics are too low to measure 
tth Yukawa coupling. At the CERN Large Hadron Collider (LHC), the process pp tth is an important search 
channel for a light Higgs mass from lOOGeV to ISOGeV [12], and the statistics will allow a measurement of 
tih Yukawa coupling at 11.9% level with nih = 120GeV [13]. Recently, at the Tevatron and the LHC, 0(as) 
QCD corrections to tth production processes have been investigated in Ref. [14, 15]. At e+e~ colliders, the 
measurement of the tth Yukawa coupling can be significantly improved through the process e+e" — > tih. With 
integrated luminosity of lOOOfb^^ , at y/s = 500GeV, the top Yukawa coupling can be measured at an accuracy of 
Sgtth/gtth — 21%, and with integrated luminosity of 500fb~^, at ^/s = SOOGeV, the measurement of top Yukawa 
coupling can reach the level of Sguh/gtth — 5% [16]. As shown in Ref. [17] that since the tth Yukawa coupling 
can be significantly different in the supersymmetric (SUSY) model from one in the SM, the measurement would 
provide a mean of discriminating between models. Thus the theoretical predictions including higher-order 
quantum effects in the different models, which also enters into the rate and dilutes the interpretation of the 
signal as the measurement of the tih coupling, should be important. At e+e" linear colliders, tih associated 
production process was discussed at tree-level many years ago [18], and its QCD corrections were also given 
in Ref. [17]. Very recently, a calculation of the SUSY-QCD corrections to the process has been presented in 
Ref. [19]. In this paper, we present the calculation of the C(QemTOj ^/m^) and 0{aemm^ i^/m^) electroweak 
corrections in the SM, the 2HDM of type II, and the MSSM. 

The paper is arranged as follow. In section 2, we present the notation and analytical calculations. In section 
3, we give our numerical results and discussions, our conclusions are given in section 4. The lengthy explicit 
expressions of the irreducible self-energies and vertices are given in the Appendix. 
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2 Evaluation 



2.1 Formalism 

Firstly, we define the kinematical variables and list a compact formula for the one-loop virtual corrections to 
e~e+ — > tth. The momenta of incoming electron e~ and positron e+ are pi and p2, respectively, and the 

momenta of outgoing top quark t, anti-top quark t and the Higgs boson h arc assigned to ki, k2 and 
correspondingly, with the center of mass energy squared s — [pi + ^2)^ = [ki + k2 + fca)^, for the process 

e"(pi) + e+{p2) ^ t{ki) + t{k2) + Kki). (1) 

Because of the smallness of electron mass, we neglect Higgs-electron-electron Yukawa coupling, the 2^3 
process (1) can be split to two parts: the electron and positron are annihilated to time-like 7 and Z, and then 
the virtual 7, Z decay to tih final states. Therefore, the differential cross section of the process (1) with averaged 
over spin of incoming states and summed over spin and color of outgoing states can be simplified, which should 
only depend on two parameters, conveniently chosen as the energy of outgoing top quark ki, and the energy 
of outgoing Higgs boson fcg, and the differential cross section including the one- loop virtual corrections can be 
written in the general form 

'^'^ ^ ^ '^'|M*'-'='=p -K2Re(M*^'='=^M+)], (2) 



dk^dk^ 4 8(27r)4 2s^^ 



spins 



where Nc is the SU{3)c color factor, M*''**** is tree level amplitude, and SM is the one-loop corrections to the 
tree-level process, which contain the irreducible corrections §][^(^^^^'^<^'^^^^'^°^) a^d the corresponding conterterms 



(self , vertex) _ 



SM = 5M""''^ -H (5M™''*^^ 4- (5M>°'' (3) 

with 



^^vertex ^ ^Mj^?''*"'^ -f- JM™''*'''^. (4) 



2.2 SUSY electroweELk corrections 

The SUSY electroweak corrections of order 0(aemTO?6/TO^) and 0(aemTOj j/m^) to process (1) arise from 
the Fcynman diagrams shown in Fig. 1. Wc carried out the calculation in the 't Hooft-Fcynman gauge and 
used dimensional reduction to control all the ultraviolate divergences in the virtual loop corrections utilizing 
the on- mass-shell renormalization scheme [21, 22]. FeynCalc [20] is used to calculate the one-loop irreducible 
diagrams. 

The relevant field renormalization constants are defined as 

tLO = {l + ^SZi)tL, tRO = {l + ^SZ'j,)tR, 

Afi J g \ 2^Z^z 1 -I- 2<5^77 J \ 

H \ f 1 + UZhh hSZHh ^ ^ ^ j 

(5) 



A\ _ f 1 + \SZaa \SZag \( a 

2^ZgA 1+2'- 



G \ UZga l + h6ZGG J \ G 
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The mass renormalization constants are defined as 



i^to = 'fn^ + Sml^, mlQ = ml+5ml, 

mto = mt + Smt, 

mlo = ml + 5ml, 

i^HO = mjj + Smjj, 

"^Ao = m\+Sm\. (6) 

As for the CP-even Higgs mixing angle a and tan/3 = V2/V1, the ratio of two vacuum expectation values of 
neutral Higgs fields, they have to be renormalized, too. We defined 

aa = {l + SZa)a, 
(tan/3)o = (l + (5Z0)tan/3, (7) 

where a and /? are not independent in the MSSM. 

Below we described in detail the renormalization constants and counterterms in the 2HDM and the MSSM, 

while ones in the SM can be obtained using standard techniques [22]. For the Higgs sector of the MSSM, using 
the results of Ref. [23] , the relevant renormalization constants are given by 

Sml = E'^'^imD-n^,, 
5ml = ^.""{mD-THH. 

d^HH = ^To Ife2=mi' 



5Z, 



hH 



QJ.2 \K^=m„, 

^ J:^\mj,)-THH 
mjj-ml 



my -ml 



5mA = T.^^{m\)-TAA, 



5Zaa 



dZcG = |fe==0, 



5Zga 



(ml) -Tag 



m^ 

m\ 

where the Higgs tadpole parameters Th, Th, Thh, Thh, Thh, Taa, Tag, Tgg, Th-h-, Th-g-, and Tq-g- 

defined as in Ref. [23] . The renormalization of CP-even Higgs mixing angle a is fixed by [24] 

dZ^ = \{5ZHh-5ZhH). (9) 

The renormalization of tan /3 is fixed by keeping the on-shell H^tVt coupling the same form to all orders in 
perturbative series [25], which can be expressed as 

5Z, = lt.nPi'-^-'-4 + '-4^-5Z,-,-+5Za-.-cotP) 

(10) 
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with 



5Z 



H-H- \k^=ml 



"Eff- 
using above the renormalization constants, for the process (1), the counterterms of Higgs self-energy can be 

expressed as 

Chh = i[k'^SZhh - rn\5Zhh - 5m\ + Thh], 
Chh = i[k^5ZHH - m?H^ZHH - Srn^j +Thh], 

ChH = i[\{5ZhH + SZHh)k'' - ^{ml5ZhH + mlSZHh) + T„h], 

Cgg = ilSZoGk"^ + Tgg\, 

Caa = i[k'^SZAA-m\dZAA- Stua + Taa], 

Cag = i[^{SZAG + SZGA)k''-^ml6ZAG+TAG]. (12) 
The counterterms of gauge boson and CP-odd scalar mixing, Z^ — G and Z^ — A are given by 

Sml 5Z^^ SZgg. 

C.G = + ^ + 

CzA = kf^m^^^Y^. (13) 
The counterterms of gauge boson and top pair interactions, ^ti and Zti are given by 

C.U = [9L{'-^+5Zi)+gi,^^-^+gU-+^^^.-^^ (14) 



where Pl,r = (1 T 75)/2, and 5e/e is absent at the order of 0{aein'm^ ^/ml,) corrections. 

The counterterms of Higgs and top pair Yukawa interactions, hti, Gti and Att are given by 

. .SZhh , SZ^ (5Zf dmt , Se Sm^l 5s^ Jsin/? Scosa 

Chu = i9uhi^ + — + ^)+9uH{— + --^- — -—^ + ^^) + 

9hJ-^]PL + - R), 

^ r L fSZcG , SZ^ , „L i^^t ^ Sci Smj SSyj 

Cgu = b«G(^ + — + — )+5«g(— + --^-— ) + 

9l^A^]PL + (i - Rl 

i 5Zaa , 5Z^ SZ^ L 6mt Se Sml Ss^ (5 cot/? 
Cau = [9nA{^ + — + —)+9uA{— + --^- — -^) + 

9^g^-^\Pl + {L^R). (15) 
Note that {i = 7, H, h, A, G, Z) appeared above are the coupling constants. 
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The counterterms of Higgs and gauge boson pair interactions, /177, h^z and hzz are given by 

^h-yz — 9zzh (-) 9^iV) 



Cu.. = [g.^Hi-^ + 5Z,,) + g^.nC-j + Igi - 2:g^ + y_ ) + (16) 
The counterterms of gauge boson and Higgs pair interactions, 7/1G, Z/iG and ZhA are given by 





5Zz^ 

= 9hGz 2 


CjhA 


5Z^_~, 


= 9hAz 2 






CzhG 


= [9hGz(. „ 



2 i5Zgg , iJ-^/j/ix , 5cw 5sw , 5s\n{(i-a) 

- H ^ 1 + ghGz{ h ■ ,a + 

^ 2 2 e c,„ sin(/5 — a) 

5ZaG , SZHh-tf, , N 

9hAz — 2 l-^ffGz^ — \{kh - kcJij., 

.5Zzz , 5Zaa , <52';i/i 5cyj Ssw 5cos{P-a) 

(-^zhA = [9hAz(.—^ 1 7, ^ -^)+9hAz(. 1 775 ^) + 

111 e Cw Sw cos(p — a) 

9hGz — ^ \-9hAz 2 iK'^h-kA)^- (17) 

With ah the counterterms fixed as above, the renormahzed amphtude for the process (1) is obtained by adding 
the counterterms to the corresponding irreducible corrections arising from the self-energy diagrams shown in 
Fig. 2-6, the vertex diagrams shown in Fig. 7 - 10, and the box diagrams shown in Fig. 11, respectively, 
which can be reduced by FeynCalc [20] and the relevant explicit expressions are given in the Appendix. Using 
above results it is straightforward to calculate the 0{aeinm'^ ^/m^) and 0{aem'm^ ^/mlj) SUSY electroweak 
corrections to process (1). 



3 Numerical results and Discussion 

In this section we present some numerical results. In the SM and the 2IIDM, we limit the lightest CP-cvcn Higgs 
boson mass rrih larger than 114GeV, while in the 2HDM of type II, the lower bound of the charged Higgs boson 
mass is 350GeV from the constraint of Bv{B X,^) [26]. In the MSSM, we use SUBHPOLE [27] to calculate 
the radiative corrections to the CP-cvcn Higgs pole masses and mixing angle a, and the CP-odd Higgs pole 
mass, which incorporates the one-loop effective potential and two-loop leading-log contributions from arbitrary 
off-diagonal stop and sbottom matrices. The inclusion of the radiative corrections in the Higgs sector is essential 
for Higgs phenomenology, because the tree-level lightest CP-even Higgs mass m;, < has been excluded by 
LEP experiments. There is ±3GeV inaccuracy in SUBHPOLE to calculate the lightest CP-even Higgs boson 
mass, due to unclear subleading-log and uncalculated higher order corrections, so we take the lower limit of m/j 
as lllGeV. We also take the lower bound of other SUSY particle masses as, m^o > 37GeV, iti^± > 94GeV, 
mj^ > 95GeV, m^^ > 91GeV [28], and Ap < 3 x IQ-^ [29]. In the MSSM, for simplicity, we assume At = Af, 
and use the relation of gaugino masses M2 = 2Mi at the electroweak scale. In all the three different models, 
the SM; the 2HDM and the MSSM, there are common features that the contribution of 7 exchange channel is 
dominant in the center of mass energy y/s region from 0.5TeV to ITeV, while the contribution of Z exchange 
channel increase with increasing -/s. 

The results in the SM are shown in Fig. 12 and Fig. 15. We present the cross section a and SUSY electroweak 
corrections to the cross sections relative to the tree-level values Acr/cr*'''"' with Act = ct^'^ — ct*''®'' as a function 
of the center of mass energy with ruh = 115GeV in Fig. 15. The corrections are at most a few percent 
and decrease with increasing ^/s. Fig. 12 gives the cross section and the corrections as a function of the Higgs 
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boson mass m/, at \/i = SOOGeV. The corrections increase as rrih increase, but still only a few percent. And we 
can sec that the corrections arc positive when y/s varies from 0.5TcV to ITcV for mh — 115GcV in the SM. 

Fig. 13 and Fig. 15 shows the results in the 2HDM. Fig. 15 gives a and Act/ct*™'' as a function of ^/s for 
tan/? = 10 and a = 0.05, assuming ruh = 115GeV, rriH = 250GeV, tua = 300GeV and mH± = 350GeV. The 
corrections imply a few percent reduction in the cross sections, and the magnitude of the corrections firstly 
increase as ^/s increase, then slowly decrease. Fig. 13a shows the dependence of a and /S.a/a^'^^^ on rrih for 
tan/3 = 10 and a = 0.05, assuming niH = 300GeV, rriA = 300GeV and mH± = 300GeV. The corrections are 
negative, and their magnitude can exceed ten percent which increase as rrih increase. In Fig. 13&, we present a 
and Acr/cT*''®® as a function of tan/3 for a = 0.05, assuming mh. = 115GeV, mjj = 250GeV, tua = 300GeV and 
mjj± = 350GeV. The corrections are about 2% reduction in the cross sections, and not sensitive to the values 
of tan (5. 

In Fig. 14 - 15, we give the results in the MSSM. Fig. 14a shows the dependence of a and Aa/a^'^°° on fj, 
for tan/3 = 44, assuming tua = 308GeV, Mi = 59GeV, At = 562GeV and mf^ = = mg = 375GeV. The 
total cross section a only slightly depend on the /i, but has a strong dependence on the cm. energy ^/s, as 
shown in Fig. 14a (left); the magnitude of the corrections Aa/a^'^^^ is a few percent in general, and can exceed 
10% for n > 700 GeV, as shown in Fig. 14a (right). Fig. 146 exhibits a and Aa/a^^'^'^ as a function of the 
lightest scalar top quark mass m^^ for tan /3 = 4, 10 and 40, respectively, with = 250GcV, /i = 220GeV, 
Ml = 65GeV, At = llOOGeV, and m^^ = Wj;^ = l.Sm,. One can see that the results are not sensitive to the 
values of tan/3, where there is a cut in the curve for tan/? = 4 with m^^ below ISOGeV due to the constraint of 
the lightest CP-even Higgs mass lower bound as shown in Fig. 145 (left). With y/s = SOOGeV, the corrections 
decrease the cross sections, and their magnitudes become smaller with mj^ increasing, but with y/s = 500GeV, 
the corrections can increase or decrease the cross sections depending on mj^, and they are less than a few 
percent in general, however, when m^^ is below llOGeV, the corrections can be over —10% with y/s = SOOGeV 
as shown in Fig. 146 (right). We show a and Acr/cr*™'' as a function of ^/s in Fig. 15 for tan/3 = 6 and 40, 
respectively, with niA = 160GeV, fj, = 220GeV, Mi = 70GeV, At = lOOOGeV, and m^-^ = m^^ = 1.5mj. The 
corrections decrease the cross sections, and their magnitudes increase as \/s increasing, which range between 
~ 8% and 4% ~ 13% for tan /3 = 6 and tan /3 = 40, respectively. 

4 Conclusion 

To summarize, we have calculated the electroweak corrections of order C(aein'Tit b/?Ti^) and 0(aem"^t b/w^) to 
process e+e" — > tih in the SM, the 2HDM and the MSSM, respectively. These corrections are a few percent in 
general, and can be over 10% in the MSSM with the lightest scalar top and bottom quark mass near the lower 
bound. The total cross sections including the electroweak corrections for process e+e~ — > tih get their maximal 
near y/s = 700GeV, and can reach 2.8 fb, 2.7 fb and 2.5 fb in the SM, the 2HDM and the MSSM, respectively. 

Note added: While preparing this manuscript three papers [31. 32, 33] appeared where the electroweak 
corrections to the same process in the SM are also calculated, the numerical results of their virtual weak 
corrections are in agreement with our results in the case of the SM. 
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V — V 



Figure 2: Gauge boson self-energy {vaVb). 



Appendix 

In this Appendix, we present the irreducible self-energy and vertex contributions. The convension, Feynman 
rules, and coupling constants agree with Ref. [34]. In the formulae followed, we have employed Passarino— Veltman 
one-loop functions Bi, Bij, {i = 0, 1) Ci, Cij, Cijk fc = 0, 1, 2) [20] and our notations agree with FeynCalc. 

The unrenormalized gauge boson— gauge boson VaVb (vaVb = w~w~ ,jj,jz, zz) self-energies as shown in 
Fig. 2 are as follow, which we only extract terms proportional to m^j,. 

^r(fc') = -^1 + 2^oK) + mlBoik^,mlml) + {m\ - ml)B^[k\mi,Ta\)\ 

y ^ = -|^[-(m?-™g)(So + 2SO + ((m?-mg)^-Hfc^K^-mg)K](fc^mg,m?) 
^r(fc') = - + 2AoK) + 2mlBo{k^ml mf) - ^m] - ml)B,{k^ mg, m?)] 

= -,^,[{Gji:;-G^:;^ + Gf:;-Gl:;^){-A,{mi) + m 

^T'^^^'^ = Y^KG,^<J"Gtr+Gf;°Gf^)(-^oK)+m2(l + B„(fc2,m2,m^)))] (18) 
with q ~ t,b. 

In order to renormalize the Higgs sector, we need to calculate Tg (s = h, H), which is the summation of the 
contribution of two tadpole diagrams in Fig. 3 (a) and (6), Tg = Ti"^ +Ts''\ where Ts°'\ Ts'^^ are expressed as 

Ti'') = -^^(Gjj + Gf;.)mfAo{m}) 

T^'^ = ^^^Gt.sAoiml) (19) 

with f = t,b and si = i^, bi [i — 1, 2). 

The unrenormalized Higgs boson— Higgs boson SaSb {saSb — hh, hH, HH, GG, GA, AA) self-energies are the 
summation of contribution shown in Fig. 4 a, 5 and c. 

5]SaSb ^ iI]''"''^(/t2)^i(S^aSb(fc2)_^j.SaSb(^2^)^5.SaSb^ (20) 
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: s = h, H 



is = h, H 



(a) 



ib) 



tadpole 

Figure 3: Tadpole diagrams with s = h, H. 




(a) 



Si 

S2 /' 

(b) 

s — s 



Sb 



Sb 



Figure 4: Higgs boson self-energy (saS&)- 



with 



"(b) 



N, 



+ (^l/;^l7i + ^/.7.^?A)(24i?oo + +6(i?i + Bn)k^)]{k^ml,m%) 



167r2 



'I67r2 



(21) 



where /1/2 = tt, bb, S1S2 = iiij,bibj, and s = U, bi = 1, 2). 




Figure 5: Gauge boson— Higgs boson mixing {vs). 
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t f t 

t - 1 

Figure 6: Top quark self-energy. 




V 



Si 



S2 



v-t-t 



Figure 7: Gauge boson— top quark— top quark vertex (vtt). 



The unrenormalized gauge boson— Higgs boson mixing {vs — zG, zA) in Fig. 5 are expressed as 



with 



N, 



+ ^hh + Gff. Gf:^ )™/. , ml , ml ) 

where /1/2 = tt, bb. 

The unrenormahzed top quark self-energy shown in Fig. 6 are expressed as 



with 



(22) 



(23) 



(24) 



(25) 



where fs = th,tH,tG,tA,bG+ ,bH+ ,x+ba,x%, = 1,2, j = 1,2,3,4). 
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The unrenormalized gauge boson— top quark— top quark vertex {vtt = •ytt, Ztt) shown in Fig. 7 (a) and (b) 



are 



Ty fefc2^PL + hl^.PL + fe7,.-PL + T^o h hl^PL + [L ^ R) 



(26) 

L 



where rf (i = 1, 2, . . . , 10) is the summation of rJ;J and rj;,^, rf = rf^'^ + rf^'^, where Tf^^^ and r^(,) 
are expressed as follow, and we get Tf^^ through the exchange of {L ^ R) in the expressions of corresponding 



(a,6)- 



2(o: 

3(a 



^ 4(a: 



^ 5(a 

pL 
6(o 

pL 

7(a: 



^ 8(a 

pL 
9(o 

pL 

10(a 

pL 
J- 1(6 

pL 
2(6 

pL 

3(6 



4(6; 

pL 

5(6; 

pZ, 
6(6 



7(b 



- 8(6 

- 9(6 



10(6 



= /2 '^f/i '^l' "^/i ^2 

4Ci2A;i-fc2 + 2(Ci+Cii)fc|)) 



-^^fif2^thGf-^ti^2 + C22) 



C2) 



= -^<^/l/2'^t7l'^/2t(^0 + Ol + C12 

^ /^Rv /^Rs r-ihs ( r~i i \ 

= I^G7.G/2''*(G/:/2-/^^o - -/.C^2 + G^;, m,,C2) 

/^Rv /^Rs r^hs ( \ \ \ 

- -^^fif2'^iAGf^t{Co + Gi + G2) 



1 



Gl,,GfpGy^'^mfiCo + 2C2) 



ig^2-«2«l-t7 

Y^G:,..G^fGj-m/(Co + 2Ci) 



/Oi; /^Rsi /-iLs2 ( 

Y^G,^,^G,-^. G^.^ (G2 



2G. 



22 J 



.^G:,,,G,^-GJ-(2Gi2 + G2) 



G:..,Gf;^G^;=fci^(Gi+2Gi2) 



167r2 "^"1 */ 



1 



G^2.xGf;^G^-f^(Gi+2Gn) 



167r2 ''^^i /* 








(27) 



with the alignments of the C function as G(A;|, (^1+^2)^, fcj, rr?^^rr?j^,rr?^^ and G(A;|, (fci+/c2)^, fc^, m^^, m^J 
for Fig. 7 (a) and (&) respectively. For both vertex ytt and Ztt, the virtual particles propogated in the loops are 
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Figure 8: Higgs boson— top quark— top quark vertex (stt). 

as follow, /1/2S ^ tth,ttH,ttG,ttA,bbG- ,bbH- , /sisa = bG-G-,bH-H-, and additfonally /1/2S = xtxtK, 
fsiS2 = xtKk,xli*Jl for 7tt vertex, /1/2S = , X?X?t;, /siS2 = thA,tAh,tHA,tAH,thG,tGh,tHG, 

tGH,x+b*J)*f^,xli*Jl for Z<t vertex, withi,j = 1,2, a, (3 ^ 1,2, k,l ^ 1,2,3,4. 

The unrenormalized Higgs boson— top quark— top quark {stt = htt, Gtt, Att) vertex shown in Fig. 8 (a) and 
(6) are 



(28) 



where vf'^ (i = 1,2,3,4) is the summation of rf^;^^ and vf^'^, vf'^ = vf^f^ + vf^'^', where Tf^^^ and Tf^^^ are 
expressed as follow, and we get through the exchange of (L ^ i?) in the expressions of corresponding 

i(a,())- 



^ l(a 



2(a 



3(a 



4(a 
1(6 



- 2(6 



- 3(b 



- 4(6 



1 



ff\'G'^^^^^i2Gj^^^mf,mf,Go + Gf^^{8Coo + 2(^2 + ^22)^? 



327r2 

4Ci2fci-fc2 + 2(Ci+Cn)fc2')) 
1 



167r2 
1 



^?G|-(G- m,,Go + G- m,,G2 + G^ ™,.G2) 



Rs 



-Gf-G-;(G- m,,Go + G- m,,Gi + G^ ™,.Gi) 



Rs 



16 



16^2^«2.i-t7 



G^!'G^,!^mfCo 



ft 



1 



167r2 */ /* 

1^^-2.1 Gfj G^-j Gi 



= 



(29) 



with the auguments of the G function as G(fc|, (^1+^2)^, fc2, m^, to^^, m^J and G(fc|, (fci+fc2)^, fc2, m'j,ml^,mlj 
for Fig. 8 (a) and (6) respectively. For vertex htt, Gtt and Att, the virtual particles propogated in the loops are 
as follow, /1/2S = tth,ttH,ttG,ttA,bbG- ,bbH- ,xtxfhl,xlx^iil, /siSa = xfKK,xli*Jl, and additionally 
fsiS2 = thh,thH,tHh,tHH,tGG,tGA,tAG,tAA,bG- G- ,bG- H- ,bH-G- ,bH- H- for /itt vertex, /siS2 = 
thA, tAh, tHA, tAH, thG, tGh, tHG, tGH for Gtt vertex, and fsiS2 t/iA, tAh, tHA, tAH, thG, tGh, tHG, tGH, 
bG-H-,bH-G- for Att vertex, with i,j = 1,2, a,^ = 1,2, k,l = 1,2,3,4, 

The unrenormalized Higgs boson— gauge boson— gauge boson {sVaVb — /177, hjZ, hZZ) vertex shown in 
Fig. 9 (a)-(e) are 



Fig^j. + T2ki^kii, + F3/c2^fc2i, + Tiki^k2„ + T^k2^ki^ + T(ie^""'^ki„k2r 



(30) 
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Figure 9: Higgs boson— gauge boson— gauge boson vertex (sVaVb). 



where Ti (i = 1, 2, . . . , 6) is the summation from ^i(a) to ri(e): = ri(a) 



+ ^He), where rj(a), ri(c) 



and rj(e) are expressed as foUow, 



1(a) 



2(a) 



3(a) 



4(a) 



5(a) 



6(a) 



1(c) 



-2(c) 



1 



72/3 



/3/1 



/2/3 



(-4C7oo - 2(^2 + C22)fci + 2(Ci + 2Ci2)fci • - 2Ciifc2) + 
^hf^^hl^tk + 0r/2G|7 )m/3(-8C„o - 2(C2 + ^22)^? + 

2(Co + Ci + 2Ci2 + C2)fci • k2 - 2(Ci + Cii)A:|) + (G^^/.^J/^GJ;;^ + Gf:fGff\Gf^}^)mf, 
(4Coo + 2C22fc? - 2(2Ci2 + C2)ki ■ + 2{Ci + Cn)kl)] 

'^^^^hh'^hh^hh + G/l/2G/3/lG/2/3^'"^lG22 + (t-/i/,G/3/,t./,/3 + Gf,f,G^-^^G^-^^Jmf,{C2 + C22)\ 

-i^i(G/i/2G/3/iG/2/3 + '^fihGf-,fGj^fJrnf,{Ci + Gn) + (G f^f^G ^-^^G ^-^^^ + G f^f^G ^-^^G ^.^^Jmf,Gn\ 
^\-yGhhGf,fGj^_f^ + Gf^fG^^^G^^j^)mf,(2Gi2 + C2) + (G f^^^G ^.^^G ^-^^^ + G f^^^G ^-^^G ^-^^.Jruf, 
(Ci + 2C12) + (Gf.^Gj- + G|- Gf- )m,3(Co + Ci + C2)] 

^^^^fihG^Gf^j^ + Gj^^G^^^G^^^^)mf,[2Gi2 + G2) + (^^^^^G^^^^G^^^^ + G/^/.G^^^^G^^^Jm/^ 
(Gi + 2G12) - (Gf, G|- G|- + G- G|- G|- )m,3(Go + Gi + G2)] 

^ \( r^Ls f^Lva r^Lvf, , ^Rs r<Rva /^^""b \^ r< j_fr<R^ f^Lva r^Lvf, y^Ls r'R-Va /^^""b \^ n 
f r-iRs f~iLva r^Rvb , f~iLs f^Rva z-tLvb \ I ^ \ r< _L M 

^ /l/2 /3/1 /2/3 /1/2 /3/1 /2/3' 13\ 



/2/3 



/1/2 /3/1 /2/3'^ 



^G^»,^G^^^3GJ^,JG2 + 2G22) 
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(a) 
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,f \ ,*■ \ 

1 - 1 

(c) (d) 
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Figure 10: Gauge boson— Higgs boson— Higgs boson vertex (vSaSb)- 



4(c 



i 5(c 
r6(c 



J- l(e 

^3(e 
^4(6 
^5(6 
^6(e 



^Z^G's^s^G^J^^Gg^s^ (Ci + 2G11) 



1 

872' 

^ ;G^s;s,GZ,,Gl,,ACo + 2Gi + 4Gi2 + 2G2) 



167r2 



GX^G:,,,i?o((fcl + fc2)2,m2^,m^J 



1 

^ 157J-2 '^S2S1 ^S1S2 






(31) 



with arguments of G functions as G(A:|, (A:i + k2)'^ , kf,m'j^,m'j^,m'j_^) and G(/c|,(fci + fc2)^, A:^, , m^^ , m^J 
for Fig. 9 (a) and (c), and we get the result of Fig. 9 (5) and (d) with the exchange of Va ^ Vb and /, si /, S2 
from the one of Fig. 9 (a) and (c). For the vertex /177, hjZ, hZZ, the virtual fields propogated in the loops 
are as follow, /1/2/3 = ttt,bbb in Fig. 9 (a), S1S2S3 = iaipi^.habpb^ in Fig. 9 (c), S1S2 = iaif3,babf3 in Fig. 9 
(e). Note that 7 couples only with scalar quarks of the same quantum numbers, i.e., jtata- 

The unrenormalized gauge boson— Higgs boson— Higgs boson {vSaSb — jhG,jhA, ZhG, ZhA) vertex shown 
in Fig. 10 {a)-{d) are 

r-a«. = ri(A:i-A:2)^ + r2(fci+fc2)^ (32) 

where ri^2 is the summation from Ti(^a) to Ti(^cL), Ti — ^i(a) +ri(fc) + • • ■ + ^i{d), where Ti(a) and rj(c) are expressed 
as follow, 

p , _ ^ r n(/^Bv i^LscL /^Rst I ^Lu r^Rsa r^Lsb \^ , ^ _ /-/^ _j_ ^ If ( rfLsa /-tLst . f^Lv /-iRsa r^Rsh \ 
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Figure 11: Box diagrams. 

+ (G^:/ + Gf^.^Gf; G|- )m,Jm,3(Co + G, + G.) - 

^^hh^hK^ll + ^hh'^hl^hlM^^ + 2G22 + 2G222)fc? + 2(12Goo + 9Gooi + 9G002) - 
6(Gii2 + 2Gi2 + G222)A:i • + 3(Gi + 2Gii + Gm + Giia)^^)] 

mf,mf,{C, + Gi - G2) - (G^:/,Gf;;^Gf;^ + Gf^^G|;;Gf;^)m/,m/3(Go - Gi + G2) ~ 

(^/i^^S.^fA + ^hh^hl^hV^-^^^^^ + 6G002 - G2fc? + 2(Gn2 - G222)fcl • fc2 + 

(Gi - Gill + Gii2)fc|)] 
-'"i(c) = ^;^Gs3SiG's2S3GsjS2(Co + Gi + G2) 

^2(c) = ^;^^s3SiG's2S3G5^j,J-Gi + G2) (33) 

with arguments of G functions as G(fc2, (fci + ^2)^, /^i , m^g, m^^, m^J and G(fc|, {ki + k2)'^ , kf , 171"^^, ml^ , mlj for 
Fig. 9 (a) and (c), and we get the result of Fig. 10 (&) and {d) with the exchange of Sa ^ Sf, and /, si /, S2 
from the one of Fig. 10 (a) and (c). For the vertex 7/1G, 7/1^, 2'/iG, ZhA, the virtual fields propogated in the 
loops are as follow, /1/2/3 = ttt, bbb, S1S2S3 = tatfft-y, babpb^. 

The analytic expressions of box diagrams in Fig. 11 are too length to present here, while we keep them in 
our numerical calculation. 
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Figure 12: cr(left) and corresponding Acr/cr*'''^'' (right) as a function of with ^/s = SOOGeV in the SM. 
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Figure 13: cr(left) and corresponding Acr/cr*'''"' (right) as a function of mh and tan/3 corresponding respectively 
to (o), and (6) in 2HDM with y/s = SOOGeV. 
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Figure 14: (T(lcft) and corresponding Acr/cr*''®'^ (right) as a function of n and the hghtest scalar top quark mass 
rrif^ corresponding respectively to (a) and (b) in MSSM with y/s = 500, SOOGeV. 
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Figure 15: C7(left) and corresponding Aa/a^""" (right) as a function of in the SM, 2HDM and MSSM. 
The other 2HDM parameters are ruh = 115GeV, niH = 250GeV, niA = SOOGeV, mH± = 350GeV, and 
a = 0.05. And the other MSSM parameters are Ma = 160GeV, At = ITeV, n = 220GeV, Mi = 70GeV and 
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